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Abstract: A new approach for continuous and non-invasive monitoring of the 
glucose concentration in human epidermis has been suggested recently. This 
method is based on photoacoustic (PA) analysis of human interstitial fluid. The 
measurement can be performed in vitro and in vivo and, therefore, may form 
the basis for a non-invasive monitoring of the blood sugar level for diabetes 
patients. It requires a windowless PA cell with an additional opening that is 
pressed onto the human skin. Since signals are weak, advantage is taken of 
acoustic resonances of the cell. Recently, a numerical approach based on the 
Finite Element (FE) Method has been successfully used for the calculation of 
the frequency response function of closed PA cells. This method has now been 
adapted to obtain the frequency response of the open cell. Despite the fact 
that loss due to sound radiation at the opening is not included, fairly good 
accordance with measurement is achieved. 
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1 Introduction 


The photoacoustic or photothermal effect is utilized for the highly sensitive 
detection of weakly absorbing sannples [U [2l |3l |4]. In photoacoustic sensors 
(PAS) radiation fronn a source, typically an infrared laser, is used to excite 
vibrational states of nnolecules. Non-radiating relaxation of these molecules 
leads to a local elevation of the temperature. Since the radiation is modulated, 
the temperature varies periodically. The temperature variation is accompanied 
by a pressure modulation, which can be detected by a microphone. As long as 
no saturation occurs, the PA signal is proportional to the concentration of the 
excited molecules. 

To obtain a large PA signal one often uses the closed sample cell as an acous¬ 
tic resonator. When the modulation frequency is tuned to an eigenfrequency of 
this resonator, acoustic modes are excited which results in signal amplification. 

Recently it has been suggested to use PAS for the non-invasive in vivo 
measurement of the glucose content in the interstitial fluid. The glucose con¬ 
centration in the interstitial fluid is correlated to the glucose level of the blood 
and the goal is to obtain a fast, cheap and painless measurement of the blood 
sugar level of diabetics [3 E]. This would also be a large step forward to a 
sufficient monitoring of blood sugar content as required by physicians. 

The T shaped resonator of the PAS for blood sugar monitoring introduced 
in [5j is depicted in Figure It consists of a broad main cylinder with a base 
which is either closed by a window or open (laser beam opening). The axis of 
the narrow resonance cylinder is perpendicular to the axis of the main cylinder. 
The microphone is mounted near it’s outer end. The probands have to close the 
absorption opening at the top of the main cylinder by pressing the skin of their 
hand to the cell. After entering the cell through the optical window respectively 
through the laser beam opening the modulated laser radiation (Uber Tuner 9, 
Daylight Solutions, California) traverses the PA cell and is directed onto the skin 
surface. The interstitial fluid which contains the glucose molecules is located 
about 50 /im to 100 /im below the surface. From here the PA sound waves are 
emitted, travel into the resonator and are detected by the microphone. 

Early tests in ambient air with such a device revealed a problem: Due to 
skin transpiration during the measurements humidity accumulated inside the 
temperature stabilized cell leading to condensation at the cell walls and the 
window. However, humidity in the sample cell increases the attenuation of 
sound waves especially at frequencies above 10 kHz [7]. One remedy could be to 
guide a gentle flow of dry nitrogen through the cell j6j. Unfortunately, this flow 
results in noise, which superimposes the PA signal. In addition, the necessity 
to amend the device by a fixture for the maintenance of a steady nitrogen flow 
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Figure 1: T shaped resonator of the PA sensor made from stainless steel. In the 
first experiments the opening at the bottom of the large cylinder (main cylinder) 
was closed with an optical window (closed resonator). The small opening at the 
top of the cell is closed by the skin of the patient. The microphone is mounted 
at the end of the narrow cylinder (resonance cylinder). 


is in conflict to the aim of a simple blood sugar sensor for domestic us. An 
alternative solution uses an open resonator and thus allows the water molecules 
to evaporate out of the cell. An open resonator has also the advantage of 
keeping the static pressure at a constant value and temperature variations due 
to skin-cell contact small. On the other hand, the opening of the resonator 
deteriorates the signal and, therefore, the detection sensitivity. 

However, a strong PA signal is of crucial importance. The determining factor 
for the signal strength is the shape of the resonator. Unfortunately, it is not 
obvious which cell geometry results in a high acoustic amplification. Testing 
different shapes experimentally is an option but would require much time and 
is expensive. It is considerably more efficient to calculate the signal strength of 
a certain geometry. A method for the analytical calculation of the PA signal 
has been available for a long time [8]. The method, however, is applicable for 
simple geometries like cylinders only. In we have combined the method of 
Kreuzer with the Finite Element Method (FEM), thus allowing to apply it to 
virtually arbitrary shapes of resonators. The procedure used to be restricted to 
closed resonators. Here, we investigate the possibility to extend the application 
to open resonators. 
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2 Theoretical background 

2.1 Calculating the PA signal for closed resonators 

The acoustic pressure p at the microphone positior0^ and a modulation fre- 
quency ou can be deternnined by solving the inhonnogeneous Heinnholtz equation 

V‘^p{r, cj) + k‘^p{f, cj) = i(jj ^ T-L{r, u) . (1) 

c is the speed of sound, k the acoustic wave nunnber and T-L{r^uj) the Fourier 
transform of the power density 7 is the ratio of the isobaric and 

isochoric heat capacities. If the absorption of the radiation by the molecules is 
not saturated and the modulation frequency is much smaller than the relaxation 
rate of the molecular transition T-L{r^uj) = al{r^uj), where X(r,cc;) is the Fourier 
transformed intensity of the electromagnetic field, a is the absorption coefficient 
for the infrared radiation in the sample. The walls of the resonator are assumed 
to be sound hard. The solution of Equation Q can be expressed as 


( 2 ) 
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where the modes^j(r) and the according eigenfrequencies cjj = ckj are obtained 
by solving the homogeneous Helmholtz equation. The modes are orthogonal and 
have to be normalized. 

The frequency dependency of the amplitudes Aj{uj) in Equation ([^ has the 
form 


Aj{ijj) = i 


AjUJ 


( 3 ) 


The contribution of a certain mode is determined by the excitation amplitude 


Aj — 


«(7 - 1) 



■XdV, 


( 4 ) 


where Vq is the volume of the PA cell and p^ the conjugate-complex pi. As 
mentioned in the introduction the sound is excited a small distance below the 
skin surface. In principle, the domain of integration has to be extended to the 
region, where the radiation is absorbed. In practice, this has been handled a 
little different (see Section 3.1). 


^In principal, one needs to integrate the acoustic pressure over the diaphragm of the mi¬ 
crophone. Nowadays diaphragms are very small and it is sufficient to calculate the pressure at 
one representative point. 
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Table 1: Used gas parameters (air). The values come from the database of the 
FE software (COMSOL Multiphysics 4, see |[9|) and correspond to a temperature 
of 20°C and a static pressure of 1013 hPa. 

density p — 1.2044 kg/m^ 

sound velocity c = 343.2 m/s 

viscosity p = 1.814 • 10“^ Pas 

coefficient of heat conduction Av = 2.58-10-2 W/mK 

specific heat capacity at constant volume cv = 7.1816 • 102 J/kgK 
specific heat capacity at constant pressure Cp = 1.0054 • 10^ J/kgK 


So far the model described above does not account for loss. If loss is small 
it is possible to incorporate the damping of sound waves through the inclusion 
of loss factors £j in the amplitude Equation Q: 


Aj{uj) Aj{uj) 


AjU 


(jj^ 


— (jj‘^ + iUJUJni 


3^3 


( 5 ) 


Various sources of loss have been identified (see, for instance HU). One ob¬ 
serves viscous and thermal dissipation in the fluid body (volume loss) and at 
the resonator boundary layers (surface loss). There is loss due to acoustic wave 
scattering at surface obstructions, due to the compliance of the chamber walls 
and dissipation at the microphone diaphragm. In particular in the region above 
10 kHz humidity inside the resonator can strongly increase the damping of 
sound waves, an effect which might be important in the context of in vivo mea¬ 
surement of blood sugar. Finally, surface roughness can increase surface loss 
due to viscosity [Jj. The combined effect of the various loss mechanisms can 
be calculated by adding the individual loss factors. 

In the present work only volume and surface loss due to viscosity and thermal 
dissipation is considered [8]|T0]. The physical parameters used here and in the 
following are compiled in Table[^ The compliance of the walls of the sample cell 
is negligible. For modern microphones, like the one used in the setup discussed 
here (Knowles SPM0404UD5), dissipation at the diaphragm is of no importance. 

Originally, the differential equation and the surface integrals appearing dur¬ 
ing the calculation of the loss factors are calculated analytically j^. Obviously, 
only simple geometries can be treated. For realistic profiles of the intensity X(r) 
the integral Q cannot be calculated without simplifications. In this paper these 
quantities are calculated by a FE software and arbitrary geometries and realistic 
intensities can be used. 
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Using the nnethod excellent agreement between numerical and experimental 
results of a T shaped test cell has been obtained [TOj. In addition, a shape 
optimization on the basis of the FE model has been performed which lead 
to the proposal of an hour glass shaped resonator with an improved signal in 
comparison to a conventional cell [T^ [T3] . 


2.2 Open resonator 

In this section the differences of an open resonator to a closed resonator are 
described. The theory comprises advanced mathematics, but we are interested 
in some basic results only m Levine and Schwinger considered the leaking of 
acoustic energy at the opening of cylindrical ducts. They assumed the sound 
wave is generated by a harmonically driven piston inside the pipe [15] ■ This 
system bears a certain similarity with our open resonator (Figure Q. If we 
ignore the resonance cylinder and regard the sound generation at the skin as 
similar to the sound generation at the piston, this becomes obvious. The results 
obtained in [14] can then be used to get an idea what is happening at the open 
end of the PA resonator. 

The i-th longitudinal eigenfrequencies of a pipe (radius i?, length L) with 
one open and one closed end is given by 


f^ = 


4L ’ 


( 6 ) 


where i = 1,2,3,.... The investigations of Levine and Schwinger revealed that 
L has to be modificated by the end correction AL = aR with a either equal to 
0.61 (unflanged case, meaning the walls of the pipe are of negligible thickness) or 
a = 0.84 (flanged case, meaning the walls of the pipe are of infinite thickness). 
The acoustic pressure drops to zero at a distance AL outside of the opening 
and the eigenfrequencies are accordingly lower as predicted by Equation Q. 
Figure [T] shows a PA cell which approximately represents the flanged case. 

The loss due to the leaking of sound energy from the laser beam opening 
can be described by an additional loss factor [7j: 


nrad _ 


1 w; 

fj Ej 


rad 


(7) 


Ej is the acoustic energy content of the resonator when the j-th mode is excited 
and is defined by a surface integral over the opening: 




op\ 


IPjf d5, 


( 8 ) 
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Re(yop) denotes the real part of the specific admittance of the opening. 

In m the low frequency behavior of the acoustic impedance (the inverse 
of the admittance) of the opening has been derived: 

zop ^ ^kR (iai + a2kR). (9) 

The wave number k is complex since the waves suffer attenuation. The numbers 
ai and a 2 depend on the frequency. In practice one often uses ai ^ 0.6 and 
02 ~ 0.25 as an approximation m- 

In leading order of kR, i.e. in the low frequency limit, the impedance is 
purely imaginary. In this limit the sound wave is completely reflected near the 
opening and suffers a phase shift of tt. When the frequency is increased sound 
energy begins to leak out of the opening. 

It seems that we could use the above equation for to obtain the admit¬ 
tance of the opening and to calculate the radiation loss factor. Unfortunately, 
the theory is applicable for u < c^crit = 27r/crit only. For the resonator consid¬ 
ered in this paper the critical frequency is about /crit = 25 kHz and much lower 
than the frequency used in experiment. Above /crit radial and azimuthal modes 
contribute in the case of a cylindrical geometry m 

The general tendency is that the amount of energy that leaks out of the 
opening of the resonator increases with frequency. However, it not possible to 
state how much acoustic power is lost in comparison to the other loss mech¬ 
anisms discussed earlier. In ^ it has been pointed out that openings not 
necessarily increase the loss substantially. 

3 Results 

In [5] the response function of the open resonator of Figure has been deter- 
mined experimentally. In addition, Pleitez and coworkers have measured the 
response of the cell after sealing the opening with a window. In both cases a 
broad band absorbing material (glassy carbon black) has been used as a refer¬ 
ence sample. The laser beam is focused to a point slightly off the sample so 
that the sample is irradiated across the whole skin opening. 

We present results for the same cells obtained from a FE model and compare 
them to the experimental response function^ As explained in the previous 
section we are not able to determine the loss factor due to sound radiation. 
Therefore, a comparison of the measured and the calculated response should 
indicate how much radiation loss in the open cell contributes. 

^The experimental data used in this paper are not identical to those used in [5]. The 
deviation is minor. 
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3.1 Closed resonator 


The FE model is well tested for the case of a closed resonator [TO] and, therefore, 
it is reasonable to begin here. All boundaries of the closed cell are assumed to 
be sound hard. 

As mentioned in Section 12.11 the source term for the excitation of sound 
vanishes everywhere outside a small region near the surface representing the skin 
of the proband or the surface of the carbon black sample, respectively. Instead 
it has been used X(r) = const in the small cylinder adjoining this surface. This 
is for technical reasons and has no significant influence on the results, as has 
been checked by varying the shape and the size of the region where the heat 
source is active. 

The microphone is mounted flush at the end of the resonance cylinder and 
the diaphragm represents its ending. Since the exact location of the diaphragm 
is not known, we have conducted a parameter study, which revealed that an 
effective length of the resonance cylinder of 7.4 mm gives the best match with 
experimental results. 

In Figure]^ the experimental response function of the closed resonator is 
depicted together with the result of the simulation. It is obvious that the 
main resonances observed in experiment also appear in the simulation. The 
resonance frequencies are predicted with very good precision (deviation < 0.2%). 
The corresponding modes are displayed in Figure For at least some of the 
resonances the width observed in experiment is larger than the width of the 
calculated resonances. This indicates loss mechanisms which are not included 
in the model (see section 2.1). 

One can speculate on the influence of surface roughness of the cell walls. 
As mentioned previously the roughness of the cell walls can increase surface loss 
due to viscosity. In cylinder cells with a length to diameter ratio significantly 
larger than one, the roughness of the surface is important for longitudinal modes 
and less important for radial modes. This is due to the fact that in longitudinal 
modes the velocity of the fluid particles is along the long cylinder barrel while 
in radial modes it is along the small basal and top area of the cylinder. In T 
shaped cells the modes cannot be classified into these categories. Still, modes 
can have a more longitudinal or a more radial character. The first and the 
third mode of Figurej^correspond to weaker resonances in experiment than one 
would expect from the calculated response. In both modes one can identify a 
longitudinal component: For the first mode in the resonance cylinder and for the 
third mode in the main cylinder. This reasoning might serve as an explanation 
for the difference in calculated and measured response functions. 

The FE model predicts a strong resonance around 50.6 kHz and a small 
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Figure 2: Response function of closed resonator. The thin line with circles show 
the experinnental result, the thicker full line the FE result. The response axis 
is in arbitrary units. The curves have been rescaled such that the resonance 
annplitudes for the resonance at about 53 kHz coincide. 



Figure 3: Dominant modes of the closed resonator. Depicted is \p\. Dark blue 
corresponds to = 0, red to the maximal value of |^|. The absolute values 
of \p\ are of no significance. The corresponding frequencies are 52.0 kHz, 
52.9 kHz, 55.4 kHz and 56.5 kHz. 
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I. 

Figure 4: Modes corresponding to resonances which do not appear in in the 
experimental response curve of the closed resonator. Depicted is \p\. The 
frequencies are 50.6 kHz and 57.9 kHz. The second mode clearly has the 
character of an azimuthal mode of the main cylinder. 



resonance at about 57.9 kHz which are not present in the measured response 
function (Figure]^. The absence of the latter can possibly be explained by a 
steep drop of the microphone sensitivity of about 20 dB in the frequency range 
55 kHz < / < 60 kHz (see Figure |^. We are not able to offer an explanation 
for the missing resonance at 50.6 kHz. Most likely some details of the cell 
geometry have not been encompassed in the FE model and this difference is 
responsible for the suppression of the resonance at the microphone. 


3.2 Open resonator 


Since we are not able to include an impedance boundary condition that describes 
radiation loss in our model we chose to model the opening in the simplest 
possible way by a sound soft boundary condition. All other boundary conditions 
and the source term are identical to section |3.1| The open boundary is not 


included in the surface integrals, which are used to calculate the surface loss 
factors [ID]. 

In Figure the experimental and the calculated response function are de¬ 
picted. Starting from the high frequency end there is a pronounced resonance 
in the simulation curve at about 58.2 kHz which is not present in the experi¬ 
mental curve. Again, this might be attributed to the microphone sensitivity as 
in the case of the closed resonator. At 52.7 and 54.6 kHz one observes strong 
resonances in the simulation result. The corresponding modes are displayed in 
Figure]^ The experimental response function shows two similar resonances at 
slightly lower frequencies (deviation about 1 %). This is not surprising since 
at the open end one has to expect an end correction which leads to larger 
wavelengths, i.e. lower frequencies. 

Without supporting the statement quantitatively it is apparent that the 
deviation of height and width between calculated and measured resonance am¬ 
plitudes is similar to the one of the closed cell. If the loss of acoustic energy due 
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Figure 5: Sketch of the response curve of the microphone (Knowles Acoustics 
SPM0404UD5). Apparently there is a substantial drop of the curve in the 
relevant frequency range from 50 and 60 kHz. 
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Figure 6: Response function of open resonator. The response axis is in arbitrary 
units. The curves have been rescaled such that the resonance annplitudes for 
the resonance at about 52 kHz coincide. 



Figure 7: Dominant modes of the open resonator. Depicted is \p\. The fre¬ 
quencies are 52.7 kHz and 54.6 kHz. 


to radiation at the laser beam opening of the resonator would be substantial, 
the accordance of the measured and calculated response should be poorer in 
the case of a open cell. Finally, the calculated response function exhibits a reso¬ 
nance at about 51.1 kHz, which has no counterpart in the experimental result. 
The discrepancy might result from the same reason, which has been offered in 
Section ItT] for the closed resonator. 

4 Conclusion 

We applied a FE model originally designed for the calculation of the microphone 
signal of PA sensors with closed resonators for the determination of the response 
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function of an open resonator. At the opening a sinnple sound soft boundary 
condition has been assunned. Connparison of nunnerical and experinnental results 
reveal a fairly good concordance. The fact that sonne nneasured resonances 
are spectrally wider than the calculated ones shows that not all relevant loss 
nnechanisnns are included in the nnodel. However, the deviation of resonance 
qualities is sinnilar for closed and open resonator. This is an indication that loss 
due to sound radiation at the open boundary is not a donninant loss effect. In 
order to get a good prediction, it is innportant to know the geonnetrical data 
with high accuracy. It nnight also be innportant to consider the roughness of 
the cell walls. Further studies should take all details of the resonator shape into 
account. 
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